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Abstract 

We compute the capture rate for Dark Matter in the Sun for models where the dominant 
interaction with nuclei is inelastic — the Dark Matter up-scatters to a nearby dark "partner" state 
with a small splitting of order a 100 keV. Such models have been shown to be compatible with 
DAMA/LIBRA data, as well as data from all other direct detection experiments. The kinematics 
of inelastic Dark Matter ensures that the dominant contribution to capture occurs from scattering 
off of iron. We give a prediction for neutrino rates for current and future neutrino telescopes based 
on the results from current direct detection experiments. Current bounds from Super-Kamiokande 
and IceCube-22 significantly constrain these models, assuming annihilations are into two-body 
Standard Model final states, such as W + W~ , tt, bb or t + t~. Annihilations into first and second 
generation quarks and leptons are generally allowed, as are annihilations into new force carriers 
which decay dominantly into e + e~ , /z + fi~ and 7r + 7r~. 
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I. INTRODUCTION 



Detection of dark matter is one of the most important - and challenging - tasks of modern 
astrophysics. A wide range of experiments have been undertaken to search for this dark 
matter directly and indirectly. Of the direct detection experiments, which hope to find the 
recoil of a Weakly Interacting Massive Particle (WIMP), the signature is the recoil of a 
nucleus after collision with a dark matter particle. The key task in these experiments is the 
distinction between the rare event (a WIMP scatter) and the common background (arising 
from natural radioactivity or cosmic rays). 

One approach is to take advantage of the motion of the Earth relative to the rotation of 
the Milky Way spiral arms. As the Earth revolves around the sun, the WIMPs generally 
scatter more frequently in underground detectors when the orbit of the Earth moves it into 
the DM "wind" that arises from the galactic rotation from the (expectedly) non-rotating dark 
matter halo. This effect raises the possibility of searching for Dark Matter through an annual 
modulation signature l|, |2j. The DAMA/LIBRA experiment has released data that displays 
an unmistakable annual modulation, consistent with this Dark Matter interpretation |3j. If 
confirmed, this would give a valuable clue to the identity of the Dark Matter that makes 
up a quarter of the critical density of the Universe. In the simplest models of Dark Matter, 
however, the DAMA/LIBRA signal is inconsistent with the the lack of observation in other 
low-background nuclear experiments, in particular XENON [4j and CDMS [5(. 

However, this conclusion does not apply to all Dark Matter models. It is possible to 
make the DAMA/LIBRA observation consistent with other direct dark matter detection 

n 

experiments if what is being observed is actually inelastic scattering |6| of the type 

X + N^ X * + N. (1) 

In this process, the nuclear recoil spectra is sensitive to the mass difference 5 = M x * — M x . 
Only WIMPs that have energies above the mass splitting can drive the up-scattering to the 
heavier state. For 5 ~ 100 keV, the recoil spectrum is sensitive to the target nucleus, raising 
the possibility that the scattering rates might differ drastically between experiments. This 
class of models has been shown to be consistent with both the recent DAMA/LIBRA data, 



as well as limits from the low-background experiments 



Matter could arise from mixed sneutrinos 
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can also be incorporated into models that explain 
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the recent tantalizing excesses from 
the PAMELA [13 1 experiment, as well as an excess from INTEGRAL via the mechanism 

n 

of Ref. [141. Should this interpretation hold, both the positive result from DAMA and the 



null results from the other low-background experiments will have been crucial in pointing us 
towards the correct theory of Dark Matter. 

Indirect detection of Dark Matter through the annihilation into neutrinos can put strong 
constraints on Dark Matter models. Determination of the neutrino rate requires knowledge 
of the Dark Matter annihilation rate in the Sun (or Earth) as well as the spectrum of 
the annihilation products. Typically, a neutrino signal is potentially observable when the 
annihilation rate is in equilibrium with the capture rate (see, e.g., 15|). This capture rate 
is determined by the interaction of the Dark Matter with nuclei. Thus, the annihilation 
rate in the Sun is very sensitive to the scattering cross section off nuclei. So, capture rates 
are tied closely to the observed rate at DAM A/LIBRA, and neutrinos from the Sun are an 
important probe of any mechanism that explains that signal. Because of the gravitational 
potential of the Sun, WIMP capture in the Sun differs in important (but calculable) ways 
from the scattering off detectors in the Earth. 

In this work, we first discuss the capture rate for inelastic Dark Matter (iDM) in the 
Sun. We then discuss the number of upward through-going muon-events that might be 
seen at detectors. This depends on the final states of the annihilation of WIMPs. The 
annihilations produce Standard Model particles, which eventually decay to neutrinos, which 
then propagate from the Sun to the detectors on Earth. A preliminary estimate of neutrino 
rates was made in [61], but here we extend the work, taking full account of the form factors 
in the scatterings, a proper treatment of the velocity distribution, as well as full propagation 
of the neutrinos from the center of the sun to the Earth. 



3 



II. INELASTIC CAPTURE RATE FOR INELASTIC DARK MATTER 



The central question in determining the rates of neutrinos is generally the WIMP capture 
rate on the sun. As WIMPs annihilate, one ultimately reaches equilibrium between capture 
and annihilation, i.e., C = 2r . Thus, the capture sets the upper bound, and often the 
expected signal of neutrinos from the sun. In this section we present analytical formulae for 
calculating the capture rate of inelastic dark matter in the Sun. 

An inelastic WIMP of mass m x can only scatter off a nucleus of mass if its energy 
E > 5(1 + m x /mjv). We will show a generalization of the results of Refs. [3] and 17], that 
allows for a discussion of the capture of inelastic dark matter. 

A. Kinematics of Inelastic WIMP scattering 

The kinematics of inelastic scattering are quite different from that of elastic scattering. 
To have any possibility of scattering at all, the WIMP must satisfy the kinetic energy 
requirement 

E mm > = s(i + H*). (2) 

The consequences of this change are purely kinematical. I.e., 

d&elastic do 'inelastic 

dEji dEfi 

but the allowed energy ranges for scattering can be quite different. In the rest frame of 
the nucleus, conservation of energy and momentum of the WIMP-nucleus scattering process 
implies a minimum and maximum amount of energy loss in such a collision: 

6 _ X 



AE mm = — + 2-^ (E - ^E(E - SX+j^j (4) 
AE max = — + 2^{E+ y/E(E - (5) 

where 

Til 

X = _ * X± = X ± 1. (6) 

m N 
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We use the parameter X to make connection with Ref. 16j . which uses fi. We use X to 
avoid confusion with the reduced mass parameter. This definition of X± is twice that for fi 
in Ref. Q. 

As a consequence, the total cross section for inelastic scattering is 



25 

^inelastic 1/1 o ^elastic- ( ' ) 

flV z 

We will typically define our cross sections in terms of 00, which is the cross section per 
nucleon in the elastic limit. 



B. WIMP velocity distribution and flux 

To compute the scattering rate of the WIMPs on the Sun, we must characterize the 
distribution of WIMP velocities incident on the Sun. We assume that the WIMPs have an 
isothermal speed distribution in their rest frame. However, this frame need not coincide with 
the rest frame of the Sun. In the Sun's rest frame, we have: 

1 ( ul 



f( u ) = ^77 7 ex P 



'o 



N(v esc )~"*-\ 

where v is the velocity dispersion, and 

u 2 Q = u 2 + v 2 Sun - 2u ■ v s , (9) 
with vsun is the speed of the Sun with respect to the WIMP rest frame. The normalization 



factor, N{v esc ) is given by Ref. [181 ]: 



N(v esc ) = (nv 2 f /2 ( erf (y m ) - ^=y m e*» ) , (10) 



where y m = v max /v, and v esc is the escape velocity. In the limit of y m — > oo this expression 
simplifies to N(oo) = (nv 2 ) 3 ^ 2 . 

We now wish to compute the flux distribution of WIMPs incident on the Sun. We begin 
by considering an imaginary surface bounding a region of radius R about the Sun, which is 



5 



FIG. 1: Coordinate system used in the derivation of the capture rate. The inset is the local 
coordinate system on the surface of the sphere, which defines the angles a and (3. 



large enough so that the gravitational field due to the Sun is negligible. For an infinitesimal 
region of area R 2 sin 9d9d<p we choose a local coordinate system (see Fig. 1) so that 

u = ucosar + usma cos j39 + wsinasin/30 (11) 
u ■ v s = uvsun(cosacos9 — sina;cos/3sin0). (12) 

Therefore the flux of WIMPs with velocity in the range (u, u + du), (a, a + da), ((3, (3 + d(3) 
through the infinitesimal area R 2 sin 6d8d<j) is 

/ u 2 + v"ij un — 2ttf g un (cos o. cos 9 — sin a. sin 9 cos /3) \ 

= exp I — J x 

\ v 2 J 

R 2 sin 9d9d(f) (—u ■ f) DM u 2 dud(cos a)df3. (13) 

If we make the approximation that v esc — > oo, then the analysis simplifies considerably. As 
we will see, we do not expect this assumption to drastically affect the results. With this 



6 



assumption is possible to do a straightforward integration over (3 (and </>) 

/ (m 2 + V 2 Sun - 2uV S un(cOSaCOs9))\ (UVSun . . Q \ 

dT = -exp I — — — — — J J ^ _ 2 sin a sin j x 

fl 2 sin 9d9 27r2nDM u 3 dud( cos 2 a). (14) 

Here, Jo is the modified Bessel function, and udm = Pdm/^x ls ^ e WIMP number density. 
We now change variables to angular momentum per unit mass = J. Therefore 



u 2 R 2 ' 

and the total flux of WIMPs becomes 



cosa= -1—, (15) 



U 2 + V 2 Sun -2uV S unCOSe^l-^ \ f Jv 

dT = exp ( — | Iq [ n „ 9 sm J x 



Rv 2 



svnOdd . DM ududJ 2 (16) 

iV(oo) v ' 

exp (_ - 2 + <n-^v Sun cose \ d(cose f* 2 n»M dudA 



J N (°°) 

where in the last line we take the R — > oo limit. Therefore as expected the flux of WIMPs 
in the forward direction is greater than that in the direction opposite to the velocity of the 
Sun. 

C. Differential capture rate 



Following the calculation in Ref. 16j, the WIMP velocity at the some r inside the Sun is 



w 2 = u 2 + v 2 (r) (18) 
and the probability for the WIMP to scatter inside a shell of thickness dr at radius r is 

n-(w)- (19) 

w 

2 ( J 2 V 1/2 
= tt-(w)-dr 1 - — - e(rw - J) (20) 
w \ r z w z J 
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Therefore the differential capture rate is a product of the flux times the probability to scatter: 
dC = ^( W )4rWp p 2 + ^ -* uv sun cos 9 ^ dcos9dr 

x 2 -^^duQ{\ mw 2 -5X + ). (21) 
iV(oo) 2 

To arrive at this expression, we performed the J 2 integral from to r 2 w 2 . The step function 
Q{\m x w 2 — <Mf + ) imposes the constraint that the energy of the WIMP is sufficient for the 
collision to occur. Therefore we find 

~dV = J x dX j Q f(^ x ) du - n vH^m x w 2 -SX + ) (22) 

where 



2irn DM ( u 2 + v 2 Sun - 2uv Sun x\ 2 



/(«, = ^P ^ ^ J u'dti (23) 

- A ™ DM ^ ■ exp f- ^H sinh ( 2 -^) u 2 du. (24) 



N(oo)uV Sun V V 2 J \ V 2 

where in the last line we have performed the integral over x. Note that if we had not taken 
the v esc — > oo limit, the upper limit of integration on the u integral would have been both a 
and f3 dependent, greatly complicating the analysis. Generalizing the result in Eqn. (A5) of 



Ref. 



16| we have 



^(w) = — ^- / d(AE)F 2 (AE)Q(AE-E 00 ), (25) 

max J AE min 

where E = mw 2 /2, E^ = m x u 2 /2, and E^fJ lc = 4fi 2 v 2 /m^ is the maximum nuclear recoil 
energy in the elastic case. F 2 (AE) is the form factor that accounts for the fact that at 
sufficiently high momentum transfer, scattering from the nucleus is no longer coherent. We 
have also defined 

_ m\m 2 N (f pZ + f n (A - Z)f 

m 2 n {m x + m N ) 2 f 2 

where m n is the nucleon mass and is the WIMP-nucleon scattering cross-section. f p>n are 
the relative proton and neutron couplings, which we take to be equal. The results for the 
differential capture rate in Eqns. (I22"l) - (l26l) represent the main analytic result of this paper. 
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FIG. 2: Contour plot of Cq/10 23 s _1 in the Mwimp — 5 plane, where vsun = 250 km/s, v = 
250 km/s, p DM = 0.3 GeV/cm 3 and <r = 10- 40 cm 2 . 

D. Inelastic dark matter capture and astrophysical uncertainties 

In this section we discuss the implications of Eqn. (1221) and the uncertainty in the capture 
rate due to the astrophysical uncertainties in vsun, v and the distribution of heavy nuclei in 
the Sun. We also consider the effect of changes in F 2 (AE) in Eqn. (125]) . The changes in 
the capture rate that result from varying these astrophysical and nuclear parameters give 
us an estimate of some of the systematic uncertainties in our calculation and an idea of the 
robustness of the current limits we discuss in Sec. IIV Al As a reference point, we use the 
values pdm = 0.3 GeV/cm 3 and o"o = 10" 40 cm 2 throughout this section: a change in these 
values will just rescale the total capture rate. Also, as reminder, this value of o"o is chosen 
because the inelastic cross-section needed to make the DAMA/LIBRA results compatible 
with direct detection experiments is of the same order of magnitude. 

Before exploring the effects of these astrophysical uncertainties, we establish a baseline 
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FIG. 3: (left) Capture rate of a 100 GeV WIMP due the different species of nuclei in the Sun, 
assuming the standard solar profile of Ref. [19|, |20(] and the same values of vs U n,v,riDM and ctq as 
in Fig. [2j (right) Expanded view of (a) for the region of S > 100 keV. 



in Fig. [2j There we show the variation of the capture rate in the my/ imp — 8 plane for the 
values Vsun — 250 km/s, and v = 250 km/s. We have also used the Helm form factor from 
Ref. [l^l for F 2 (AE). As in the case of elastic dark matter, the capture rate decreases with 
increasing WIMP mass. This is in part due to a simple decrease in the number density of 
WIMPs. However, there are two additional effects: it is more difficult for heavy WIMPs 
to lose energy in collisions with the relatively light nuclei in the Sun, and as the mass of 
the WIMP increases, it becomes more difficult to satisfy the minimum scattering energy 
condition of Eqn. (T5|). The minimum scattering energy condition in Eqn. (T5]) leads to a 
suppression of the capture rate due to lighter nuclei, as shown in Fig. [3j 

To focus on the effect that the inelasticity has on WIMP capture, in Fig. [3] we show the 
capture rate of a 100 GeV WIMP on each species of nuclei in the Sun. We assume the 
same values of vs U n,v,nDM and <7o as in Fig. El We see that for 5 > 30 keV, the scattering 
of WIMPs off hydrogen and helium is highly suppressed, and in the range 100 keV < 5 < 
150 keV the capture due to scattering off iron dominates the other elements by a factor of 4. 
As a reminder, this range of values 100 keV < 5 < 150 keV provides the best fit to the current 
DAMA data. The importance of the heavy elements for the capture rate is not unique to 
the inelastic case: for both elastic and inelastic (spin- independent) models scattering off the 
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FIG. 4: The suppression of the capture rate in the inelastic case relative to the elastic scattering 
case with comparable gq shown as a function of Mwimp for two choices of S = 100, 150 keV. 

heaviest elements can dominate. (This is not true of spin- dependent scattering.) This is 
perhaps counter-intuitive, as the abundances of the heavy elements are substantially less 
than hydrogen and helium. However, these tiny abundances are compensated by the fact 
that the (coherent) spin-independent couplings scale as the square of the atomic number. In 
fact, when one accounts for additional kinematic factors, for large Dark Matter masses the 
spin-independent cross-section scales as the fourth power of the atomic mass number. These 
A A factors can more than overcome the relative scarcity of heavy elements. Some reviews in 
the literature dangerously neglect the contributions of these heavy elements for simplicity. 

Depending on the choice of 5, there is actually a slight enhancement of capture due to 
iron relative to the elastic case. The reason is that form factor suppression is smaller for the 
inelastic case compared to the elastic one: for the same loss in WIMP energy AE the recoil 
energy of the nucleus is Er = AE — 5 for the inelastic scenario compared to En = AE for 
WIMPs that scatter elastically. The net effect of inelastic scattering compared to the elastic 
case is shown in Fig. HI 

We now turn to sources of uncertainty in the capture rate. The importance of heavy nuclei 
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FIG. 5: (a) Standard Solar profile from Ref. Qjl^l- (b) Constant profile, (c) Linearly decreasing 
profile, (d) Effect of changing the profile of the heavy elements in the Sun on the capture rate of 
inelastic WIMPs for S = 100 keV with the remaining inputs the same as in Fig. 

in the capture of inelastic WIMPs implies a sensitivity to the abundance and distribution 
of these nuclei in the Sun, which is somewhat uncertain. To quantify this dependence on 
leavy nuclei in the Sun we consider two variations to the Standard Solar 
201 ]: a uniform distribution as a function of radius for all heavy nuclei 
and a linearly decreasing distribution with the same average number density as the uniform 
one (which we view as somewhat of an extreme case). We show these distributions, along 
with the resulting capture rates in Fig. [51 Again, we have set 5 = 100 keV and values of 
vsun,v,nDM and do that are the same as those in Fig. [21 From Fig. [5] we see that if the 
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FIG. 6: Dependence of the total capture rate for 5 = 150 keV on vsun and v. The red (solid) curves 
bound the region swept out by Cq for the range of 200 km/s < vsun < 300 km/s and the blue 
(dashed) band corresponds to the values of Cq for the range of 200 km/s < v < 300 km/s. 

number density of heavy nuclei in the center of the Sun is increased the overall capture rate 
can be significantly increased. Also, the mass fraction of iron in the Sun will directly affect 
the capture rate. 

The threshold for scattering in Eqn. ([2]) also suggests that the capture rate is sensitive 
to variations in the astrophysical inputs vs U n and v, as well as perhaps the escape velocity. 
First, we address the issue of the finite escape velocity. We do not expect the capture rate 
to be very sensitive to this quantity. This is because the particles that are most easily 
captured are not those particles in the Boltzmann tail. These particles have the most energy 
and are harder to capture. This is in strict contrast to scattering at DAMA/LIBRA, where 
only the most energetic particles in the Dark Matter distribution are able to drive the 
inelastic transition, simply because these particles do not have the benefit of mining the 
gravitiational potential energy of the Sun. (The intuition that the result is insensitive to v esc 
was crudely confirmed by varying the (large) limit of the u integral when performing the 
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numerical integration, and seeing that the capture rate was insensitive.) In Fig. Owe show 
the predicted range of values 200 km/s < vs U n < 300 km/s with a central value of 250 km/s 
in the red (solid) band and 200 km/s < v < 300 km/s with a central value of 250 km/s 
in the blue (dashed) band, where we have used the Standard Solar Model heavy nuclei 
number densities, Pdm = 0.3 GeV/cm 3 and a = 10~ 40 cm . If one assumes no co-rotation of 
the WIMP halo, then the velocity of the Sun relative to the WIMP halo is determined as 
Vrot + ^peculiar with \v pecu i iar \ ~ 20 km/s [2l|. Recent measurements place v rot = 250 km/s. 
We take this as a central value, but allow a range values to account for the possibility of 
errors in its determination and some co/counter- rotation. If one further assumes that the 
halo is virialized then v = v rot , but we allow v and vs U n to vary independently in the plot. 

The variation in the capture rate illustrated in Fig. [6] is mainly due to the fact that mostly 
only slow moving WIMPs can lose sufficient energy to be captured. This effect dominates 
that of the WIMP minimum energy scattering condition in Eqn. ([2]). Therefore increasing 



v has the effect o 
Refs. 



decreasing the capture rate. This result is the same as that observed in 
16| and [17|]. Even varying over this generous range, only a factor of two change in 
the capture rate is observed. 

Finally, we note the importance of form factors to the discussion at hand. Due to the 
energy threshold in Eqn. (J2j) the WIMPs undergoing capture are relatively energetic and 
have to transfer a substantial portion of their energy to the nuclei. Therefore the recoil 
energy of the nuclei is large, and the process is very sensitive to the form factors used. If 
the Helm form factor overestimates the suppression for iodine at DAMA relative to iron for 
capture, it is possible that a relative factor of 0(1) could result. 

All together the uncertainties in astrophysical and nuclear inputs could lead to a factor 
of a few in the uncertainty of the capture rate of WIMPs in the Sun. The DAMA/LIBRA 
modulation signal is much more sensitive to the velocity distribution of the Dark Matter 
particles, on the other hand. This is because scatterings there only arise from particles in 
the tail of the velocity distribution. It is encouraging that the results for capture presented 
here are so robust. 
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III. CAPTURE VS. ANNIHILATION IN THE SUN 



We now review the interplay between capture and annihilation in the Sun. When the 
capture and annihilation rates are sufficiently large then equilibrium will be reached between 
the two processes. 

Unless otherwise stated, we will assume that the capture and annihilation rates are 
sufficiently large so that the WIMPs are in equilibrium, and thus the annihilation rate 
is just one-half the capture rate: Ta = C/2. 

A. Thermalization in Elastic Models 

A potentially major difference between the inelastic picture and the conventional picture 
is the question of what happens to WIMPs once they have scattered the first time and been 
trapped in the Sun. For conventional WIMPs, repeated scatters cause the WIMP to settle 
into thermal equilibrium in the center of the Sun. In contrast, for an inelastic WIMP, if the 
kinetic energy is suitably low, it is possible that no subsequent inelastic scatterings can take 
place. We must consider what occurs under these circumstances as well. 

Let us begin by assuming that there are only inelastic scatterings. In this case, the WIMP 
will proceed through the Sun, scattering off iron (which is abundant at the 10~ 3 level), until 
such a point as it has insufficient kinetic energy to scatter off iron. There will be a range 
of orbits of varying ellipticity. Some will be circular, but most will be elliptical, and thus 
proceed more deeply into the Sun during their orbits. To determine how extensive this set 
of orbits will be, we must consider how much kinetic energy a given WIMP will have in the 
interior of the Sun. 

A WIMP starting at rest from the surface of the Sun will have a velocity of approximately 
1240 km/s in the center of the Sun. This is more than enough kinetic energy to scatter off 
iron for the parameters of interest. In fact, using the density profile in the solar models of 

nn . 

|19L 1201] . simply going from R = 0.2Rq to the center will give enough energy to inelastically 
scatter off iron. Thus, most WIMPs will be contained within this region. 

We can now estimate whether the density in this region will be high enough to bring the 
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system into equilibrium (i.e., the outgoing annihilation rate equals the capture rate). If we 
assume a capture rate Cq, then (neglecting annihilation) the total number of particles in 
the Sun is Cqt , where r is the age of the Sun. If the WIMPs are captured within a radius 
R = ei?Q, the present annihilation rate is then 

r ann « (27) 

The annihilation will reach equilibrium with the capture if 

™ * " 5 x 10 m r {— ) y • < 28) 

The thermal cross section needed to produce the correct relic abundance (3 x 10~ 26 cm 2 ) is 
sufficent to put the WIMPs into equilbrium unless it decreases with velocity (e.g. is p-wave). 
Should the WIMPs fail to reach equilibrium, there is an approximate overall suppression of 
the annihilation rate given by the LHS divided by the RHS. Note that even in models where a 
p-wave cross section determines the relic abundance, there is also a subdominant s-wave cross 
section. This s-wave scattering could well be sufficient to put the WIMPs in equilibrium. 

Up to this point, we have considered particles with only inelastic scatterings, but (sub- 
dominant) elastic scatterings are common in many explicit models of inelastic Dark Matter. 
Tor instance, mixed sneutrino and neutrino models have sizable Higgs mediated couplings 
(| (a n > 10 -45 cm 2 ). Higgsinos-like models have loop-mediated contributions which can be 
important 22j (a n > 10~ 48 cm 2 ). Generically, inelastic models will have cross sections which 
are o n ~ {$/ Mw imp) 2 ® inelastic ~ 1CT 52 cm 2 [l0|. So how large a cross section is needed to 
thermalize? 

The total number of scatters in the Sun over the age of the Sun for a WIMP would be 
P ~ i w ( P \ ( ° 



N scat * m a n VTKlx ^ 1500kg . m - 3 J ^xlO-^cm 2 ^ 1,300 km- iT 1 J (29) 
where p is the typical density seen by the WIMP and m p is the proton mass. The solar density 
in the core is typically 1.5 x 10 5 kg • m~ 3 , roughly 100 times larger than the mean solar density. 
To thermalize via scatterings off of protons, the WIMP must scatter of order Mx / m p times. 
Thus, a trapped 500 GeV WIMP scattering in the interior (with high densities) would need 
a cross section of approximately 10~ 48 cm 2 to thermalize. 
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Models with light mediators typically will not have this, but they often have the 
Sommerfeld enhancement to boost annihilation rates to equilibrium. Models that are not 
pure SU(2) doublets, but instead have some mixing with SU(2) singlets after EWSB generally 
have large enough Higgs couplings. Finally, models which are pure SU(2) doublets - without 
Higgs couplings - generally have large enough loop induced cross sections to thermalize them 
in the core, although this is marginal for WIMPs much heavier than 500 GeV. 

Thus, we believe even in the inelastic case, it is reasonable to generally consider models 
that annihilate in equilibrium with their capture rates. We will make this assumption in 
what follows. 



IV. NEUTRINO RATES 



Armed with a capture rate, we can now predict the high energy neutrino flux (and the 
resulting muon rates) on the Earth. These rates will depend on the products of the WIMP 
annihilations. We considered the final states of WW, ZZ, bb , ti, rf, cc, and light jets. We 
will not explicitly show separate results for WW and ZZ as they are very similar. 

If WIMP annihilations are solely to charged leptons of the first two generations the 
annihilation products interact strongly enough in the Sun that they typically come to rest 
before producing neutrinos. The neutrinos from such decays are low energy, and will not 
give an observable signal. A similar statement is true for particles that annihilate to first two 
generation/gluon jets. The neutrino production comes from fragmentation to heavy quarks, 
and is very suppressed. For our numerical results, we follow the work of Ref. 23], which 
builds on the work of Ref. [241 ]. 



A. Existing Limits 

Over the region of interest, the strongest bounds are placed by the Super-Kamiokande 
experiment [25] and the recent data from IceCube-22 [26]. At present, the limits are 
calculated assuming specific annihilation channels, which complicates the extraction of limits 
for a new model. 
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To place limits on the model at hand using the data from Super-Kamiokande, we calculate 
the number of expected signal events at the detector, following the methods of Ref. 27]. We 
first compute the flux of neutrinos from the Sun based on the appropriate annihilation 
channel. 

We then propagate the neutrinos from the center of the Sun to the Earth, using the 



formula of Ref. 



231 ] which are available at Ref. 



These results for propagation are 



consistent with those of Ref. 
Super-Kamiokande by using the formula 

'da 



291 ] . From this neutrino flux, we calculate the rate of muons at 



dE^dE„A eff (E^ 



Nevts = T 

The muon effective area A. 



dE, 



-P P + 



da u 



dE, 



-Pi 



d<t> 
dE, 



R[i{E^ 



(30) 



e ff(E^) = 1200 m 2 , and the livetime, r, is given by | x 1670 days 
(the one-half is to account for night-time). The muon range, R{E P ) ) is approximately given 
by 

R»(E») = -^log a+ i E " , (31) 



— log — , 

p/3 a + (3E thresh 

where p is the relevant density, a ~ 2.0 MeV cm 2 g _1 . (3 varies depending on the material. 
For standard rock, (3 ~ 3 x 10~ 6 cm 2 g _1 , whereas for water, (3 ~ 4.2 x 10~ 6 cm 2 g -1 . At 
Super-K, there can be conversions both in the nearby rock and in the water. For simplicity, 
we will use the value 3 x 10 -6 cm 2 g" 1 , as rock usually dominates. If we instead use the 
value for water the event rate increases by roughly 30%. To compute the number events in 
Eqn. (1301) we use the following neutrino-proton and neutrino-neutron cross-sections [301 ]: 



da 



up 



dE„ 



da,- 



dE„ 



TV 

2 

7T 



-a < /. y dpu ~~r~ Op U ^2 



-m n G 



m p G 2 F 



, h K 

+ nv El 



(32) 
(33) 



Here a pu = 0.15, b pu = 0.04, a nu 
anti-neutrinos can be found by a. 



•25, b nu - 
= b b - 

u nvi u pv 



0.06, and the corresponding expressions for 



bpy i b nu 



a, 



pv 



The results from Eqn. (1301) can then be compared to the data for upward going muons 



coming from the Sun as reported in 25j. To calculate the number of observed (N b s ) and 
expected background (N^ d ) events as a function of the WIMP mass we use data on upward 
going muons from Fig. 5 of Ref. 25j. The region of the sky surrounding the sun used in the 
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search varies as a function of WIMP mass. We use Fig. 8 of Ref. 31] to specify the angle 
about the sun that Sun as a function of the dark matter mass, through going muons is to 
Wimp annihilations in the Sun clS 8b function 

Kt = NT + OMg* (34) 
where the factor of 0.9 takes into account the fact that only 90 % of the signal is contained 



with the window angle in Fig. 8 of Ref. 3lJ. Therefore the maximum allowed number of 
events is the value N^ 9 ^ al needed so that the cumulative poisson distribution function is 10%, 
for N b s number of observed events and number of expected events. The breaks in the 
curve around 90 and 225 GeV correspond to places where the size of the cone around the sun 
changes. In that cones size, there is a fluctuation in the number of events observed, which 
affects the limit. Note that the optimal cone size should actually depend on annihlation 
channel, as the neutrino spectra (and hence the correlation with direction) change with 



final state. However, the cone size in 3jJ is calculated instead with a "breadbasket" final 



states[4l|], so there is some uncertainty on the exact position of the curves. 

If the Dark Matter annihilates to a hard channel (e.g., W bosons, top quarks or tau 
leptons) then the limits from Super-K are very constraining. For cross sections consistent 
with the DAMA result, the typical event rates at Super-K would be too large by some two 
orders of magnitude (see Fig. [7]). If, instead the Dark Matter annihilates to a softer channel 
(bottom quarks or charm quarks), then the tension is lessened. 

Nevertheless, it is fair to say that even in the case where the Dark Matter annihilates 
through a relatively soft channel such as charm (or more so with bottom quarks), there is 
tension with the existing limits from Super-Kamiokande. At larger masses, the inelastic 



explanation for DAMA comes into tension with results from the CDMS experiment 



32|. 



However, N-body simulations have found significant structure in the high-velocity tails of 



velocity distributions 331 ]. and these may open the ranges of parameter space significantly 



when included [341 ] . 



It is important to note bounds from some channels (e.g. tau leptons) are much stronger 
than those from others (e.g. b quarks). The result is that the dominant annihilation mode 
may not provide the most stringent bound on a given Dark Matter candidate. The plots in 
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FIG. 7: We plot the maximum allowed cross section per nucleoli (in units of of 10 -40 cm 2 ) times 
the branching ratio to the given annihilation channel. The curves are derived by imposing that the 
number of expected events (background + expected signal) is consistent with the number observed 
at the 90% CL. For reference, cross sections consistent with the DAMA result are typically of the 
size a > 2 x 10 -40 cm 2 , with 100 keV < 5 < 140 keV . We consider annihilations to several two- 
body Standard Model final states: W bosons, top quarks, charm quarks, b-quarks, light quarks, 
and tau leptons. For large masses > 250 GeV in the hard channel Ice-Cube22 stronger bounds (see 
text and Fig. El). 



Fig. [7] can be used to extract the maximum branching ratio to a given channel. 

Note, however, that at higher WIMP masses, (and especially for hard annihilation 
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FIG. 8: As in the previous Figure, we plot the maximum allowed cross section per nucleon (in units 
of of 10 -40 cm 2 ) times the branching ratio to a given annihilation channel. The curves are derived 
by imposing that the number of expected events (background + expected signal) is consistent with 
the number observed at the 90% CL. We also show limits extracted from the recent IceCube-22 
results, which are (subtantially) stronger for the hard channel in the regime where they apply. 

channels) the strongest limits come from the recent results of the 22-string run at IceCube. 



In Ref. j26|, the collaboration quotes the maximum allowed values of the solar annihilation 
rate (= one-half capture rate in equilibrium) as a function of the neutralino mass for two 
different annihilation modes: XX — > bb (soft) and XX — > W + W~ (hard). For the WW 
annihilation mode, IceCube-22 excludes capture rates above Cq > 1.2 x 10 22 s _1 (3.2 x 10 21 
s _1 ) for m = 500(250) GeV. No limits are quoted below 250 GeV. Comparing these rates 
with the capture rates of the previous section, we can see that these limits are very strong. 
For the softer annihilation channel, the limits are weaker: Cq < 2.8 x 10 23 s _1 for m = 500 
GeV. No constraints on lower masses for soft channels are given. For reference, at 500 (250) 
GeV, a typical capture rate (for a = 10~ 40 cm 2 ) is 8(20) x 10 23 s _1 . Thus for b quarks, the 
IceCube constraint on the cross-section is roughly comparable to that of Super-Kamiokande 
for masses greater than 500 GeV (see Fig. [8]). Because of the large energy threshold at 
IceCube-22, it is dangerous to extrapolate to lower masses, or to draw a strong conclusion 
about annihilation to charm quarks. It is likely that IceCube should be able to probe these 
modes soon for higher masses (see next section). 

For annihilation to W bosons, IceCube-22 gives a bound that is nearly a factor of 4 
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stronger than Super-Kamiokande for masses greater than 250 GeV. Again, it is dangerous 
to extrapolate to lower masses. Something similar would be expected from other hard 
annihilation channels (tau leptons, top quarks, Z bosons). It is clear that IceCube-22 places 
strong constraints on the hard annihilation channels. 



35H to derive a 



B. Future Telescopes 

We now turn to the sensitivity of future neutrino detection experiments to neutrino 
signals. We concentrate on the case of the IceCube experiment. We calculate both 
the expected number of events induced by neutrinos in the detectors, and the expected 
backgrounds that arise from the flux of atmospheric neutrinos. 

For the atmospheric neutrino backgrounds, we use the data from Ref. 
power law fit, which takes the approximate form 

<$> atm (E u ) « 5.8 x KT 2 ^ 314 cm~ 2 GeV" 1 sr" 1 sec" 1 

§ atm (E v ) « 5.6 x 10- 2 E^ 3 - 20 cm" 2 GeV" 1 sr" 1 sec" 1 (35) 

An experiment can focus only on the direction of the Sun to significantly reduce this 
background. We integrate over a 3° region for IceCube. 

To calculate the actual number of events at IceCube, we again use Eqn. fl30l . For 
the purposes of determining the energy loss in ice, we take a ~ 2.0 MeV cm 2 g _1 and 
(3 ~ 4.2 x 10 -6 cm 2 g -1 . We take a threshold energy of -E^res/i = 50 GeV, which is 
somewhat optimistic but may be attainable, especially with the planned inclusion of the 
of the DeepCore array j^i]]. We take the effective area A e ff(E lx ) for IceCube from Ref. {37 ]. 
We find approximately 40 background events arising form atmospheric neutrinos within a 3° 
window. The corresponding number of signal events, normalized to a cross section o"o = 10 -40 
cm 2 are shown in Fig. [9j 

If nature has chosen the past of inelastic Dark Matter, perhaps the mostly likely reason 
that it has gone undetected thus far is annihilation proceeds to modes that provide soft 
neutrinos. Thus, it will be important for IceCube (and other future experiments such as 
ANTARES) to try and push their energy threshold as low as possible. 
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FIG. 9: The number of signal events expected in one year of running at IceCube, assuming a cross 
section of do = 10 -40 cm 2 per nucleon and a 5 = 125 keV and annihilation to W bosons (left) and 
b quarks (right). Note such a large cross section is excluded in the case of the W boson, see Figs. [7] 
andU 

V. CONCLUSIONS 

Inelastic dark matter provides an exciting proposal to reconcile the results of 
DAM A/LIBRA with other direct detection experiments. However, its large cross section 
makes detection through neutrino telescopes a particularly important constraint. We have 
found that these experiments can place strong limits on these models, and so it is important 
to clearly state what those constraints are. 

Under the assumption that inelastic dark matter s-wave annihilates with a thermal 
cross section, it certainly appears that for hard annihilation channels (e.g. tau leptons, 
W bosons, monochromatic neutrinos) such a scenario is excluded. For other cases where 
the annihilation products are softer sources of neutrinos (charm quarks or b quarks), the 
answer is less clear. Annihilation to charm, as well as other light quarks seems safe, at 
least within astrophysical uncertainties (though the tension in the charm channel is large at 
larger masses). Annihilation into b quarks seems excluded at the factor of two level (and 
more for higher masses). When including the variety of uncertainties, both from astrophysics 
as well as other issues such as nuclear form factors this is borderline. It is important to note, 
however, that annihilations into Higgs bosons generally include a component of r as well, 
which even as subdominant contributions are often the dominant limit. In particular when 
considering the IceCube-22 limits, particles above ~ 250 GeV are strongly constrained. 
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We should note that these limits are possibly evaded in elastic models when the 
annihilation rate is p-wave suppressed, and equilibrium between annihilation and capture 
has not yet been achieved. 

Importantly, models in which the dark matter annihilates into new, light force carriers that 
dominantly decay into e + e~, fi + n~ and 7r + 7r~ appear safe from these constraints, because 
muons and pions stop before decaying and producing neutrinos. 

Another possibility to evade these bounds is that the Dark Matter might dominantly 
annihilate to large multiplicity final states. In this case, the neutrino energies are degraded, 
as the energy is shared amongst a large number of final decay products, and limits might 
easily be evaded. This could occur naturally in light mediator models, where one might have 
XX — > 00 — > 4b, or models where the Higgs boson dominantly decays via pseudo-scalars jssj], 
in which case XX —>■ hh — > 4a — > 86. If the light states (0 or a) are allowed to decay to r 
leptons or directly to neutrinos, then tension may still exist. A detailed examination of such 
decays and related model building is left for future work[39j. Because the bounds on the 
channels with energetic neutrinos are so much stronger, it is possible (or even likely) that 
a sub-dominant decay mode with hard neutrinos could provide the strongest constraint for 
models of this type. This is also of potential relevance if Higgs boson decays are involved 
and both r leptons and 6-quarks are potentially present. 

In the case where Dark Matter is not its own anti-particle, and possesses a conserved 
quantum number, the Dark Matter abundance is due to a small excess of Dark Matter 
over anti-Dark Matter. Then captured Dark Matter may not annihilate and there will be 
no signal. Such a scenario is incompatible with potential annihilation signals observed by 
PAMELA and FGST. 

Ultimately, the space of models to explain DAMA through inelastic scattering is still 
large, but the space is strongly constrained by these neutrino telescopes. Should future 
direct detection bear out the presence of inelastic WIMPs, particularly at higher masses, 
these null results should allow us to distinguish among a variety of candidates. 
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